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Perturbation theory for non-abelian gauge theories at finite temperature is plagued by infrared 
divergences which are caused by magnetic soft modes ~ g 2 T, corresponding to gluon fields of 
a 3d Yang-Mills theory. While the divergences can be regulated by a dynamically generated 
magnetic mass on that scale, the gauge coupling drops out of the effective expansion parameter 
requiring summation of all loop orders for the calculation of observables. Some gauge invariant 
possibilities to implement such infrared-safe resummations are reviewed. We use a scheme based 
on the non-linear sigma model to estimate some of the contributions ~ g 6 of the soft magnetic 
modes to the QCD pressure through two loops. The NLO contribution amounts to ~ 10% of the 
LO, suggestive of a reasonable convergence of the series. 
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Figure 1: Feynman diagram contributing to the pressure. 

1. Introduction 

Perturbation theory for static quantities of non-abelian gauge theories at finite temperature 
(T) features three scales: %T associated with the non-zero Matsubara modes that arise due to the 
compactness of the Euclidean time direction, and gT,g 2 T, where g is the gauge coupling, which 
are associated with the screening of colour-electric and magnetic fields, respectively. A prohibitive 
obstacle for the evaluation of perturbative series are the well-known infrared divergences due to 
the magnetic modes. For example, the thermodynamic pressure corresponds to all vacuum loop 
diagrams evaluated with finite temperature Feynman rules, and one (Z + l)-loop example is shown 
in Fig. [jj The Matsubara sum also contains a bosonic zero mode, and its contribution to the pressure 
is, parametrically, 

~g 2 ' (t J d 3 p^j + p 2, (p 2 + m 2 )- 31 ~ g 6 T\g 2 T /m) l -Hov I > 3, (1.1) 

where we have introduced a mass as infrared regulator in the propagators. Clearly, for m — > 
the contribution diverges. One may argue that the full theory dynamically generates a magnetic 
mass scale, m = const. g 2 T , to regulate this divergence. After all, we know that non-perturbatively 
the theory is finite. However, in that case we see that the coupling constant drops out of the 
expansion parameter and all higher loop orders contribute to the order g 6 in the pressure. This is the 
Linde-problem [1], which occurs sooner or later in any perturbatively computed observable, e.g. in 
the Debye mass already at NLO, no matter how weak the coupling. Thus, at finite temperature, 
perturbation theory is only well-defined to some observable-dependent low order. Note that this 
IR-problem in the zero-mode sector is equivalent to the one in 3d Yang-Mills theory. 

In this contribution we reconsider a resummation method for 3d Yang-Mills theory, viz. the 
finite T zero mode sector, which has been developed a while ago [2,3]. In view of current and future 
heavy ion collision experiments there is an urgent need for non-perturbative theoretical predictions. 
On the other hand, straightforward Monte Carlo simulations of lattice QCD do not work for finite 
baryon densities or dynamical problems involving real time (for a recent review, see e.g. [4]), 
motivating also analytical attempts. As the simplest observable to test our resummation scheme, 
we therefore consider the contribution of the infrared sector represented by the 3d Yang-Mills 
theory to the thermodynamic pressure of the QCD plasma. Of course, the equation of state can 
be computed without major problems on the lattice, at least for a temperature range T <5T C . The 
idea here is to try and develop resummation schemes, test them against the pressure and in case of 
satisfactory results apply them to other quantities. In the thermal context, resummation schemes 
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have already been successfully employed in scalar theories through four-loop level [5], but it is 
not straightforward to generalise these methods to gauge theories. Recently a scheme based on 
Hard Thermal Loop terms was attempted for non-abelian gauge theory, though not including the 
magnetic modes [6]. We therefore briefly review the contribution of the 3d gauge theory to 4d 
thermodynamics, before presenting our resummation scheme and evaluating it through two-loop 
order. 



2. High T effective theory from dimensional reduction 

At sufficiently high temperatures, the previously mentioned scales get separated hierarchi- 
cally, g 2 T <C gT <C 7lT. This suggests to successively integrate out those scales, which has been 
done systematically at the two-loop level. Here we follow [7], to which we refer for details and 
references. 

In the first step, the hard non-zero Matsubara modes are integrated out, leading to so-called 
EQCD as an effective theory for modes ~ gT and softer, 



Pqcd(T) = PE (T) + ^\n J DA a k DA a exp(- J d d x££ E 



= ^TrF 2 + Tr[D k Ao} 2 + m 2 JrA 2 +X E l \TrA 2 ) ) 2 + X E 2) TrAt ) + ..., 



where to leading order the parameters of the effective theory are p E ~ T, m E ~ g T z ; g E 
g 2 T; X E ^ ~ g 4 T; X E ^ ~ g 4 T. Note that by performing the Matsubara sum the theory has 
become effectively three-dimensional, and the electric gauge field component now represents an 
adjoint scalar field. The dots represent higher dimension operators that are suppressed by inverse 
powers of T. This effective theory still contains the two dynamically generated scales gT and g 2 T. 
In a second step, the Aq field and all modes living on the scale gT can be integrated out as well. 
This step already requires some resummation, but can still be performed as a series in the coupling, 

j DA^DAq exp (— J d d x££ E ) = p M {T) + ^\n J DA a k exp(- J d d x££ M ) 

J£ M = l -TrF 2 + ..., (2.1) 

with matching coefficients pm ~ m E T and g 2 M ~ g\. This is an effective theory for the ultra-soft 
modes ~ g 2 T, its leading term being equivalent to 3d Yang -Mills theory. 

At this stage we can evaluate the contributions to the pressure from the different momen- 
tum scales separately and with appropriate methods. The contributions Pe(T) +Pm(T) can be 
computed as a series in the coupling constant, whereas the contribution of Pm(T) is completely 
non-perturbative. This is the part of the theory giving the g 6 contribution to the pressure where the 
Linde-problem surfaces, and for dimensional reasons one can give its dependence on the coupling 
constant, 



p G (T) = £ln| DA a k exp(-S M ) ~ Tg 6 



Pg(T) 

— a r$ 8m 



Tlx -2e dACA (Anf r u \e ' ""'2m G 



a G [ I + 81n-^-)+j3 G + 0(£) 



(2.2) 
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The coefficients (Xg,Pg first receive contributions at the four-loop level. There are similar (and 
more lengthy) expressions for Pe(T) + pm(T), for the full result see [7]. It is well-known that 
whatever renormalisation is necessary and sufficient in the vacuum will also be sufficient at finite 
temperature. The epsilon poles appearing in the contributions of the different scales are due to the 
break up of the momentum integration range, thus the sum of all divergences has to cancel in the 
full pressure. 

Here our interest is merely in the contribution from the 3d Yang-Mills part. Whereas CCg is 
known analytically [7], the infrared coefficient j8g is receiving contributions from all loop levels 
and needs to be evaluated non-perturbatively, e.g. by lattice simulations of the effective theory [8], 
leading to a numerical value (with a certain error bar). In the following we discuss possibilities for 
an analytical evaluation by resummation methods. 



3. Resummation schemes: general idea 

The general idea of a resummation is, at some given order in a perturbative scheme, to sum 
up higher order contributions (infinitely many in our case) into the current one. In order to avoid 
double counting, these contributions must then be left out at the order where they naturally occur. 
In other words, the perturbative scheme gets reorganised in some systematic way. (For a discussion 
of various schemes along those lines used in the context of thermal field theory, see [9]). This can 
be formalised by rewriting the Lagrangian serving for the perturbative expansion as [10] 

1 



Sf M (V£K)+^(V£X)-£Sf^V£X) 



(3.1) 



where X stands for a generic field; J£m is the theory under study, here taken as the 3d Yang-Mills 



theory of Eq. (2.1); the modification J^, contains fields of the original Lagrangian and possibly 
auxiliary fields; and I is a counting parameter introduced to systematise the resummed expansion, 
which is now an expansion in powers of I. At the end of a calculation it is to be set to i = 1, for 
which the Lagrangian is identical to the original one and the theory is unchanged. However, during 
a perturbative evaluation the contributions of Jz?0 are subtracted out one order higher than they are 
added in, hence at any finite order of the expansion the result will differ from the unresummed one. 
For example, the one-loop vacuum diagrams contributing to the pressure will be of order ft while 
two-loop contributions count as ft. In particular, if J^, is chosen to represent a mass term for the 
gluon, this will regulate the infrared divergences. 

Simply adding and subtracting a mass term for the gluon will not be sufficient, however, as 
we would obtain gauge-dependent and hence arbitrary results at every finite order in our resummed 
scheme. The question then is how to resum in a gauge-invariant way. The answer is that we need to 
resum such that we maintain the Slavnov-Taylor identities of our gauge theory. Clearly, this is not 
just to make things pretty, but a necessity in order not to change the original theory in the process 
of resumming. 



4. The non-linear sigma model providing a gluon mass term 

A well-known and time honoured way to introduce mass for gauge fields is by means of the 
Higgs effect, i.e. gauged sigma models [11]. For example, in the SI/(2)-Higgs model used in the 
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electroweak theory, we have 

j&f ff = Tr(D,<I>) t (D i <I>) + i u 2 Tr(<I> t <I>)+2ATr(<I> t <I>) 2 (4.1) 

= \( + iT a K a ) with T a = tr 
Y 2 V ' 2 

o 2 = v 2 + (n a ) 2 with <0>= v . 

The gauge field gets a mass m = gMv/2 through the expectation value of the scalar field. 

However, in the current context we don't want an additional physical scalar particle, so we 
decouple the a by sending }i , X —> °° while keeping v constant, which takes us to the non-linear 
sigma model. Choosing a parametrisation by unitary matrices generalises to SU (N) and we have 

S>t=Tr[(D4)*D4], <Hx) = — exp^T^r), 7i' a - % a + 0(n\n 5 , ...), (4.2) 

gM m 

with the usual covariant derivative D, = <?; — igM^i. What is the meaning of the remaining scalar 
degrees of freedom? This can be understood by considering the resummed partition function [10] 



Z = J DAD0 A FP exp 



- j (Sym + S ( j ) —£S ( j ) +S g f) 



(4.3) 



with some gauge fixing that only depends on the gauge field, and the corresponding Faddeev-Popov 
determinant App. We can now transform to unitary gauge by A, — > Af with U = exp(in a T a ) and 
integrate over the auxiliary field. Up to a delta-function this gives 



DQA FP e- s *f = 1, (4.4) 



such that we end up with 

Z = jDAexp -j(S Y M-m 2 J Tr A 2 + £m 2 J Tr A 



(4.5) 



which is just pure gauge theory with a gauge invariantly resummed mass term. Unitary gauge 
is not well suited for perturbative higher order calculations because of the bad UV-behaviour of 
gauge field propagators. Hence we shall proceed in a slightly different way later on, but these 
manipulations illustrate that the auxiliary field is merely a gauge degree of freedom and does not 
add anything undesired to the theory. 

Obviously, the choice of the resummation term above is not unique. There are other 
choices compatible with gauge invariance that have been considered in the past, like = F^-^jF^, 
with F^ = £u a BF a R [10]. This term does not require an auxiliary field, but at the price of being 
non-local. Another non-local choice is based on the Chern-Simons eikonal and appears naturally 
in the framework of Hard Thermal Loop resummations [12]. Yet another possibility to give mass 
to the gauge field in a gauge invariant way is by the pinch technique [13]. Different gauge invariant 
additions/subtractions correspond to different ways of resumming the theory, and it is difficult to 
say a priori which ones are better or worse. Moreover, there is no small expansion parameter here. 
Writing m = Cg 2 M = Cg 2 T, the dimensionless expansion parameter (modulo factors of 1/ (47r)) is 



2 2 

Sm _ 8 m 



m c Sm 



(4.6) 
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ref. 


m /§M 


1-loop gap eq. 


[12] 


0.38 




[2, 10] 


0.28 




[13] 


0.25 


2-loop gap eq. 


[14] 


0.34 


lattice Landau 


[15] 


0.456(6) 



Table 1: Comparison of magnetic mass values for the pole mass of transverse gluons in SU(2), m = Cgj^ 
as obtained from gap equations and from gauge fixed lattice simulations. 



where C is a numerical coefficient that will be discussed in the next section (see also Table [p). The 
coupling constant drops out and we are effectively expanding in a dynamically generated number. 
Thus, there is no limit of parameters (not even gM — > 0) in which the scheme is guaranteed to work. 
The convergence properties are due to the dynamics of the theory and can only be seen empirically 
by computing several orders of a given quantity. 



5. Self-consistency and the magnetic mass 

A priori the coefficient C for the infrared cut-off is not determined. Since it gets subtracted 
out again at higher orders in the resummation scheme it could, e.g., be treated as a variational 
parameter with some freedom to optimise the expansion. On the other hand, the resummation 
scheme of choice can be used to compute any observable, in particular one might consider the self- 
energy of the gluon. Defining the "magnetic mass" to be the pole mass of the transverse part of the 
full gluon propagator, we obtain a gap equation for m, 

^ 2) = ? +r #-n^w~J^ for p2 = ~ m2 

n trans (p 2 = -m 2 ) (l + = -m 2 )) = 

Note that the pole of the transverse self energy is gauge-invariant order by order, i.e. this definition 
of m is gauge invariant. Results for the solution of the gap equation using various resummation 
schemes and gauge fixed lattice propagators are summarised in Table p]. Expectedly there is some 
scatter in the values, but it is not much more than the two-loop collection computed with the scheme 
based on the non-linear sigma model [14], which amounts to about 15%. This suggests that the 
schemes might possibly lead to a reasonable convergence. 

There is also an evaluation of the magnetic mass in a gauge invariant lattice calculation, based 
on the following observable [16] 

i(w w \TjAd ( \ / 77 w sT exp[-(3m-2m)|x-.y|] , (5.1) 

( ( F ij) ( x ) U ab (*>?) \ F lm) GO ) 

where U^(x,y) is a Wilson line in the adjoint representation. In perturbation theory, the correlator 
in the numerator is at large distances dominated by three-gluon exchange, whereas the denominator 
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Figure 2: Lowest masses from the exponential fall-off of the field strength correlators in Eq. (5.1). Both 
masses are diverging at the same rate, but the difference can be extrapolated to the continuum. 



is dominated by two-gluon exchange, hence the ratio should fall off with one gluonic pole mass. 
In a non-perturbative reasoning these correspond to gluelump correlators, i.e. glue with certain 
quantum numbers bound to a static adjoint source. The ratio of these correlators measures the 
smallest excitation energy with the quantum numbers of a gluon. Note that, due to the diverging 
self energy of the adjoint Wilson lines, neither gluelump correlator has a continuum limit, but the 
ratio does and the energy difference is finite. The resulting continuum extrapolation is shown in 
Fig. g, and gives m = 0.36(2)g^. 

Finally, the mass parameter in the gluon propagator is not a physical observable and not the 
goal of our investigation, but merely a quantity that appears as a regulator for our resummation 
scheme. The test of the latter is in computing a physical observable. In the context of thermal field 
theory, this has been done by using the linear sigma model as a dimensionally reduced version of 
the electroweak theory in order to calculate the electroweak phase transition as a function of the 
Higgs mass. Standard perturbation theory to leading order predicts a first order phase transition 
for all Higgs masses, and cannot be extended beyond LO because of the massless W-bosons in 
the symmetric phase. Using the resummation similar to the one described here, the critical Higgs 
mass for which the electroweak phase transition disappears and turns into a smooth crossover was 
predicted [2, 17] and is within 10% of the corresponding lattice results [18, 19]. 



6. Application to the pressure 

Let us now apply the screened perturbation theory based on the non-linear sigma-model to the 



calculation of the pressure of the soft magnetic field modes, Eq. (2.2). In the framework of the 
resummation outlined above, the Lagrangian in the exponent of the first line is supplemented by 
the sigma-model term, the counter term as well as gauge fixing and ghost terms. How many powers 
of the auxiliary fields K a {x) are needed depends on the order of the calculation. At the two loop 
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level, or to order £ , a four-point vertex is the highest that is needed, and we have 

+ i (d,-^) 2 + ^m 2 ifif + (a«(c fl )*)^c a + ^V'c" 
+ 5M^r foc AfA5^A^| + ^glf£f abe f cde AlA b jA'jA'j 

+ \gM^£f abc {diK a )A b K c + g M V£f abc (d i (c a r)A b c c - l -g M ^mf bc K a c b {cy 

+ -2m1 (2 5 ab 5 cd + d abe d cde\ ^'^c^ b^ d 

8 m 2 \N J 

-^g 2 M £% (J_8«b 8 cd + d abe d cde _ f abe f cde^ ^ay^b^c^ 

- X -m 1 £A a i A a i - X -^m 2 l% a % a - t,m 2 £c a *c a . (6.1) 

Note the additional Feynman rules that result from this effective Lagrangian. Besides the vertices 
for the auxiliary fields and their interactions with gluons, ghosts and themselves, there are also 
two-point vertices associated with the counter terms, cf. Fig. ^| (left), since they are treated as 
interactions. Through order £ , i.e. two loops, we then have to compute the diagrams in Fig. || 
(right). 



p 
p 

a .. X ... b •< > V£{f) = im 2 eS ah . 

p 

Figure 3: Left: New counterterms induced by the effective Lagrangian Eq. ( |6.l[ ). Right: Feynman diagrams 
for the contribution of magnetic QCD to the pressure at order £. 

Our calculation is performed in a general -gauge. This is an excellent way of checking the 
calculation - the pressure as a physical quantity has to be independent of the gauge parameter £ . 
For intermediate checks, it is also useful to observe that the set of diagrams has three obviously 
gauge invariant subsets. The one-loop diagrams from the first line, and the two-loop diagrams just 
correspond to an unresummed non-linear sigma model which must be gauge invariant order by 
order. Correspondingly, the counter term diagrams in the last line are also a gauge invariant subset. 
According to the reshuffling of the expansion, they appear only at order £ although they technically 
are one loop diagrams. 



o * o ♦ ; 
CO + © + 



i ) 



O * O + O 
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The results for these three sets of diagrams for SU (N) are (using g 2 M = g 2 T , cf. Section ^ 



PGi0 (T)=C\N 2 -l)^T 4 
PgAT) = -£C 3 (N 2 -l)^T 4 
Pg > £> = £ C 2 N(N 2 -\)-^- 



21/1 jl \ 3 21 21, 2 

-64U +41n 2C^J + 16-32-T6 ln 3 +0(£) 



(.6.2) 



Because of the resummation, we are now getting a contribution of order g 6 already to leading order 
in this scheme, whereas in unresummed perturbation theory this would occur at four loops. We 
then extract the following coefficients (with numerical values for £ = 1, N = 2), 

a G = -^^k 2 £ = -1.015582 

1 0S /~|3 

Pg/> = — ^ 3 = 3 - 630132 
c 3 

Pc,ct = "64^ 7T 3 = -5.445198 

Pg# = -(y + 21 In ^^tt 2 = 0.196301 (6.3) 

In order to evaluate these contributions, a number for C had to be specified. To be fully self- 
consistent, we have taken this to be the magnetic mass evaluated in the same resummation scheme 
for the case of SU(2), i.e. C « 0.28 [2]. We now have two consecutive orders for j3 and may 
get a first glimpse of the convergence properties. While the resummed perturbation is organised 
in orders of £, a convergence check by comparing different orders in I would not seem to make 
much sense. By construction, the counter term diagrams get subtracted one order in £ higher 
than the other diagrams with the same number of loops and the same integral structure. Hence, 
when asking for the convergence properties, it would seem natural to proceed loopwise. Thus, in 



Eq. ( |6.3[ ) the contribution of all one loop diagrams is /3 G p + fic,ct = — 1.81, while the genuine two- 
loop contribution is j3 G £i , which is about 10% correction. We would expect the counter terms with 
two loops, entering at Z 2 -level, to be of the same order of magnitude as the two-loop contributions 
evaluated here. Thus, comparing the one-loop and two-loop results appears to be promising in 
terms of convergence. We are presently performing the three-loop computation in order to check 
this behaviour. 

Another interesting observation is the fact that if we take C to be the coefficient of the magnetic 
mass, evaluated in the same resummation scheme, then C ~ N. Thus, ,/V drops out of the expressions 



in Eq. (5.3) and j8g becomes Af-independent, as expected on general grounds [7]. 



7. Conclusions 

We discussed the possibility to gauge-invariantly resum the 3d Yang-Mills theory such as 
to self-consistently include a dynamically generated pole mass for the gluons, which serves as a 
regulator for infrared divergences. In such a scheme a perturbative expansion to arbitrary orders 
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is possible. However, the dimensionless expansion parameter is a dynamically generated number, 
hence there is no parametric limit in which convergence of the series is guaranteed. Rather, the 
convergence properties can only be inspected after the calculation of several orders. We have 
applied a resummation scheme based on the non-linear sigma-model to evaluate the contribution 
of the soft magnetic gluons to the pressure of the QCD plasma. Comparison between one- and 
two-loop contributions suggests that there is hope for a reasonable convergence of the series. 
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